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Abstract 



We report an X-ray study of the evolved Galactic supernova remnant (SNR) G156.2+5.7 based on six 
. pointing observations with Suzaku. The remnant's large extent (100' in diameter) allows us to investigate 

its radial structure in the northwestern and eastern directions from the apparent center. The X-ray 
spectra were well fit with a two-component non-equilibrium ionization model representing the swept-up 
^ ' interstellar medium (ISM) and the metal-rich ejecta. We found prominent central concentrations of Si, S 

I I and Fe from the ejecta component; the lighter elements of O, Ne and Mg were distributed more uniformly. 

2 ' The temperature of the ISM component suggests a slow shock (610-960 km s~^), hence the remnant's age is 

estimated to be 7,000-15,000 yr, assuming its distance to be ^1.1 kpc. G156.2-I-5.7 has also been thought to 
. emit hard, non-thermal X-rays, despite being considerably older than any other such remnant. In response 

to a recent discovery of a background cluster of galaxies (2XMM J045637. 2+522411), we carefully excluded 
its contribution, and reexamined the origin of the hard X-ray emission. We found that the residual hard 
X-ray emission is consistent with the expected level of the cosmic X-ray background. Thus, no robust 
evidence for the non-thermal emission was obtained from G156.2-f5.7. These results are consistent with 
' the picture of an evolved SNR. 

Key words: ISM: abundances — ISM: individual (G156.2+5.7) — supernova remnants — X-rays: 
ISM 



1. Introduction lists G156.2-I-5.7 as a mixed-morphology (MM) SNR for 

this reason, it does not fit the definition of MM SNRs in 

G156.2+5.7 (RX04591+5147) was initially discovered several ways: the X-ray shell is easily distinguishable and 

in X-ray with ROSAT (Pfcffermann ct al. 1991). Its age there are neither CO clouds nor OH (1720 MHz) mascrs 

was estimated at ^26000 yr, based on the Scdov model, in its vicinity. A recent optical observation suggests that 

^ . hence is thought to be an evolved supernova remnant the distance to G156. 2-1-5. 7 is smaller than that previously 

(SNR). The distance is comparatively nearby (~1.1 kpc; estimated (Gerardy & Fesen 2007). Their main argument 

Pfeffermann et al. 1991) and its large apparent size (^100' is that several bright Ha filaments correlate with nearby 

in diameter) enables us to study the detailed plasma struc- interstellar clouds (~0.3 kpc). However, the X-ray and 

ture. While G156.2+5.7 is one of the brightest SNRs in the radio observations both agree that G156.2-I-5.7 is lo- 

X-rays, its radio emission is quite faint; this is its most cated at 1-3 kpc, i.e., an evolved SNR (Pfeffermann et al. 

remarkable feature. Shortly after its discovery in X-rays, 1991; Reich et al. 1992). The distance and the age of this 

the radio continuum emission was also detected with the SNR remain open questions. 

Effclsberg 100-m telescope (Reich et al. 1992). As they The X-ray spectrum of G156.2+5.7 is also largely typ- 

point out, the surface brightness at 1 GHz (5.8x10"^'^ ical of an evolved SNR. Yamauchi et al. (1993) observed 

W m~^ Hz"-'^ sr"-'^) is the lowest among all the known G156.2-I-5.7 with the Ginga LAC (large area proportional 

SNRs. The radio morphology along the northwest and counters) and detected extended X-ray emission in the 1.2- 

southeast rim shows limb brightening (Reich et al. 1992; 3.0 kcV energy band. Yamauchi et al. (1999) performed 

Xu et al. 2007) that is typically seen in so-called "barrel- an additional observation with ASCA and found that its 

shaped" SNRs (Kesteven & Caswell 1987). In contrast spectrum consists of two components with different fea- 

to its radio morphology, G156.2-I-5.7 has centrally-filled tures. The soft component has several emission lines, 

thermal emission in X-rays. While Lazendic et al. (2006) suggesting a thin thermal origin whereas the hard compo- 
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Fig. 1. Left: ROSAT PSPC image of G156.2+5.7 with the Galactic coordinate overlaid. The FoV of Suzaku XIS are shown with 
solid squares. Middle: Radio (92 cm) image of G156.2+5.7 obtained by the Westerbork Northern Sky Survey (WENSS) with the 
Westerbork Synthesis Radio Telescope (WSRT). The coordinate scales are matched in both images. Right: True-color XIS image of 
G156.2-I-5.7. The energy bands of red, green and blue correspond to 0.5-0.7 keV, 0.7-1.0 keV and 1.0-5.0 keV, respectively. The red 
circular regions were excluded from our spectral analyses. 



Table 1. Summary of the 6 Suzaku observations 



Name 


Obs. ID 


Obs. Date 


R.A., Dec. (J2000) 




Effective Exposure 


NW_rim 


501075010 


2007-Feb-16 


04'^56'"54".9, 52°24'12" 


.6 


50.5 ks 


Center 


501106010 


2007-Feb-17 


04^58^54^.9, 51°44'36" 


.1 


51.2 ks 


E_rini 


501074010 


2007-Feb-18 


05'^03'"13".5, 51°37'06" 


.2 


53.3 ks 


East 


504082010 


2010-Feb-21 


05h01"W.O, 51°41'41" 


.6 


50.3 ks 


NW 


504081010 


2010-Mar-03 


04h57"M2M, 52°04'31" 


.1 


52.9 ks 


S_rini 


504080010 


2010-Mar-04 


04'^59"M1".8, 51°00'36" 


.0 


52.6 ks 



nent is featureless. However, due to the lack of statistics, 
they were not able to decide whether the featureless hard 
component has a thermal origin or a non-thermal origin. 
Katsuda et al. (2009) observed two rims (northwest and 
east) and the center with Suzaku. They found two ther- 
mal components, one of which they ascribed to interstel- 
lar material heated by the forward shock, and the other 
reverse-shock heated ejecta. They also detected hard X- 
ray observed from the northwestern rim and the center, 
which they concluded is likely to be non-thermal syn- 
chrotron emission. If this is the case, G156.2-t-5.7 would 
be the oldest SNR producing non-thermal X-ray emission. 
However, a recent Suzaku observation discovered a cluster 
of galaxies (2XMM JG45637.2-H522411; Yamauchi et al. 
2010) located in the northwestern rim. Since the surface 
brightness of this extended source is higher than that of 
the rim, we need careful consideration of its contribution 
to determine the true origin of the hard component. 

Katsuda et al. (2009) also estimated the progenitor 
mass to be less than 15Mq from the metal abundances 
of the thermal emission detected in the central por- 
tion. While no compact source has been discovered in 
G156.2-I-5.7 (see Kaplan et al. 2006), their result suggests 
that its origin is a core-collapse supernova (SN) rather 
than a type la SN. More recently, Hudaverdi et al. (2010) 
observed the center of this remnant with XMM-Newton 
and measured its metallicity. From the results, they sug- 



gested that G156.2+5.7 is an 0-rich SNR. 

We performed 6 pointing observations of G156.2+5.7 
with Suzaku. The fields of view (FoV) cover three sec- 
tions of the rim and three inner regions, and include the 
ones discussed in Katsuda et al. (2009). In this paper, 
we report the ejecta distribution and review the origin of 
the hard X-ray emission detected by Katsuda et al. (2009) 
along with the newly observed regions. We also discuss 
the age of this SNR based on these analyses. 

2. Observations and Data Reduction 

We summarize the 6 observations in table 1. Data were 
taken in 2007 and 2010, during the Suzaku AOl and A04 
phases, respectively. The observed regions are shown in 
figure 1. The squares represent the FoV of the Suzaku 
XIS (X-ray Imaging Spectrometer; Koyama et al. 2007). 
The data were analyzed with version 6.9 of the HEAsoft 
tools. For data reduction we used version 9 of the Suzaku 
Software. The calibration database (CALDB) used was 
the one updated in March 2010. We used revision 2.2 
of the cleaned event data and combined the 3x3 and 5x5 
event files. In order to exclude background flare events, we 
obtained the good time intervals (GTIs) by including only 
times at which the count rates in the various detectors are 
less than +3cr of the mean count rates. 
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Table 2. Surface brightnesses of hard emission in G156.2+5.6 and two nearby observations 



Name 


I 


b 


surface brightness at 2.0-10.0 keV 




[deg] 


[deg] 


[10^-^^ erg cm^^s^-^deg^^] 


G156.2+5.7 NW.rim 


155.5 


5.8 


1.88+0.26 


G156.2+5.7 E.rini 


156.7 


6.1 


1.61+0.23 


G156.2+5.7 S_rim 


156.8 


5.3 


1.56+0.22 


IRAS 05262+4432 BGD 


165.0 


5.6 


1.78+0.25 


SGR 0501+4516 BGD 


161.4 


1.9 


1.89+0.26 




Radius [arcmin] 



1-n fluctuation of tile CXB SB 



CXB SB [10-" ergs cm-^ s ' deg^] 



Fig. 2. Left: Surface brightness (SB) distribution of the CXB component. Filled small circle, open circle, cross, and square 
correspond to the results at E_rim, East, Center, and NW, respectively. Green (left) corresponds to the results at S_rim. Black and 
Red triangles (right) correspond to the results in the NW_rim, except in the regions of 1'.2 and 4'.0 around 2XMM J045637. 2+522411, 
respectively. The vertical error bars correspond to Icr errors. Right: Histogram of the surface brightnesses. The horizontal arrow 
represents the 1 a spatial fluctuation of the CXB surface brightness estimated by Tawa (2008). 



3. Spectral Analysis 

For the foUowing analysis, the energy range of 0.4-8.0 
keV was used for XISl (back-illuminated CCD; BI CCD) 
and XISO and XIS3 (front-illuminated CCDs; FI CCDs), 
since there was no signal above 8.0 keV after the back- 
ground subtraction (described below). In order to gener- 
ate response matrix files (RMFs) and ancillary response 
files (ARFs), we employed xisrmfgen and xissimarfgen 
(Ishisaki et al. 2007). We used XSPEC version 12.6.0 
(Arnaud 1996) for all the following spectral analyses. 
We excluded 2XMM J045637.2+522411 and three other 
point-like sources identified from the XIS image (0.5-5.0 
keV). The extracted regions are shown in figure 1. 

3.1. Background Subtraction 

Since G156.2+5.7 is a large diffuse source and the FoV 
is filled with the SNR's emission, it is hard to obtain 
the background spectra from inside the FoV. We also 
have no background data from the neighborhood of the 
G156.2+5.7. For the following analysis, we used a gen- 
erated non X-ray background (NXB) spectra by employ- 
ing xisnxbgen (Tawa et al. 2008). We also took into ac- 
count the contributions of the local hot bubble (LHB) 
and the cosmic X-ray background (CXB) by assuming 
literature-based models. For the LHB, we used an APEC 
model (Smith et al. 2001) with solar abundances (Anders 



& Grevesse 1989) and an electron temperature of 0.1 
keV (Miller et al. 2008). For the CXB we employed a 
broken power-law model with photon indices of r=2.0 
(<0.7 keV) and r=1.4 (>0.7 keV) as shown by Miller 
et al. (2008). Hereafter, we call this model the CXB 
model. We also reviewed the effect of the Galactic Ridge 
X-ray emission (GRXE). The GRXE surface brightness 
at I = 100°, |6| < 1° is ~ 1.4 x IQ-^^ erg cm-^s-^deg-^ 
in the 0.7-10.0 keV band (Sugizaki et al. 2001). Since 
G156.2+5.7 is located outside the FoV of Sugizaki et al. 
(2001) and away from the galactic plane, the contribution 
of the GRXE is expected to be negligible. For exam- 
ple, GRXE surface brightness above 2.0 keV decreases 
about one-fortieth from 6 = 0° to 6 ~ 9° (Revnivtsev 
et al. 2006). Additionally, the average 0.7-10.0 keV 
band surface brightness of G156.2+5.7 is ~ 4.3 x 10~^ erg 
cm~^s~^deg~^ (Yamauchi et al. 1999) which is an order 
of magnitude higher than the minimum Suzaku surface 
brightness of ^ 3.4 x 10~^° erg cm^^s^^dcg^^ that we see 
in the S_rim observation. In both cases, we conclude that 
the effect of the GRXE is negligible in G156.2+5.7. 

Table 2 shows the surface brightness of the CXB and the 
GRXE near G156.2+5.7 and within a few degrees using 
observations of IRAS 05262+4432 (narrow-line Seyfert 1; 
Kawaguchi et al. 2008) and SGR 0501+4516 (soft gamma 
repeater; Enoto et al. 2009). Since both are point-like 
sources, we excluded the source regions with radii 4'.0 
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and fitted the spectra obtained from the remaining regions 
with the CXB model. Comparing the surface brightness 
of each observation, the level of the CXB and the GRXE 
is uniform within the range of error. Figure 2 left shows 
the surface brightness distribution of the CXB component 
inferred from the fit to each annular region in G 156. 2 +5. 7 
(the details will be argued hereinafter). The filled cir- 
cles, open circles, crosses, and squares correspond to re- 
gions from the E_rim, East, Center, and NW, respectively. 
Green represents regions from the S_rim. The result sug- 
gests that the level of the CXB (and the GRXE) inside 
this remnant is also unchanged within the range of error. 

3.2. Model Fit for the Rim Regions' Spectra 

Since the rims are mostly dominated by an interstellar 
medium (ISM) component (Katsuda et al. 2009), we ap- 
plied a single-component VNEI (NEI vcr.2.0; Borkowski 
et al. 2001) model with TBabs (Tiibingcn-Boulder ISM 
absorption model; Wilms et al. 2000) in XSPEC for the 
rim spectra. In this model, the abundances relative to 
the solar values (Anders & Grevesse 1989) of N, O, 
Ne, Mg, Si, S and Fe were allowed to vary freely and 
that of Ni were linked to Fe. The other element abun- 
dances were fixed to their solar value. Other parameters 
were all free; the electron temperature kTg, the ionization 
timcscale r (a product of the electron density and the 
elapsed time after the shock heating), and the emission 
measure (= / rieriYidl, where and rin are the number 
densities of hydrogen and electrons and dl is the plasma 
depth). We also allowed the column density A^h to vemis- 
sion measureary freely. For the hard X-ray tail Katsuda 
et al. (2009) applied the CXB model plus an additional 
power law component. They used the same value of A^h 
for the CXB as for other components. In our analysis, it 
is fixed to 3.6xl02icm~2 calculated by nh FTOOLi and 
separated from the A^h for G156.2-I-5.7. The additional 
power law component was not applied for our analysis 
unless the CXB surface brightness exceeded the literature 
value. 

Katsuda ct al. (2009) found excess emission over the 
CXB above 3 kcV in the NW_rim and argued that the ob- 
served excess originates from a non-thermal emission from 
G156.2-I-5.7. In their analysis, they excluded a region of 
1'.2 in radius where 2XMM J045637.2-H522411 is located. 
Although Yamauchi et al. (2010) estimated the exten- 
sion of the X-ray emission from 2XMM J045637.2-h522411 
to be at least ^50" (the surface brightness within it is 
-2.64 X 10-1° erg cm-^g-^deg-^ in the 2.0-10.0 keV band) 
based on the XMM-Newton data, we noticed that its ex- 
tent might be larger. Figure 3 left and right panels show 
a full-band (0.3-8.0 keV) and a hard-band (2.0-8.0 kev) 
image, respectively. In figure 3, we show the small cir- 
cular region of 1'.2 in radius as used by Katsuda et al. 
(2009). We find that the X-ray emission from the clus- 
ter extends beyond this circle even in the hard band. 
Figure 4 shows a radial profile from the center of 2XMM 
J045637.2-H522411 in the 2.0-8.0 kcV band. The dotted 

^ http:/ /heasarc. gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh. pi 




Fig. 3. Full-band (0.3-8.0 keV; left), and hard-band (2.0-8.0 
keV; right) images of NW_rim. The FoV of the XIS are shown 
with white squares. The radii of small and large circles are 
1'.2 and 4'.0, respectively. 




2 4 6 8 10 

Radius [aicmin] 



Fig. 4. Radial profile of NW_rim in the 2.0-8.0 keV energy 
band. The center corresponds to the position of 2XMM 
J045637. 2-1-522411. The dotted vertical line represents the 
radius of 1'.2. 

vertical line corresponds to the radius of 1'.2. Figure 4 
shows that the surface brightness decreases to half at the 
radius of ^^3'. To rule out the possible contamination from 
2XMM J045637.2-h522411, we therefore excluded a circu- 
lar region with radius 4'.0 as shown with the white large 
circle in figure 3. The fitting result after excluding this 
larger region is shown in figure 5. The best- fit parameters 
are shown in table 3. 

The surface brightness of the hard component in the 
NW_rim is calculated to be 1.88x10^^^ erg cm^^s^^dcg^^ 
(2.0-10.0 keV; table 3). On the other hand, the CXB 
surface brightness is estimated to be 1.57±0.32xl0-ii erg 
cm-^s-^deg-^ based on blank sky observations with the 
XIS (Tawa 2008). In their analysis, a spatial fiuctuation 
14% (Icr) is considered for the XIS FoV (Tawa et al. 2008). 
Accordingly, the surface brightness of the hard component 
in the NW_rim (and also E_rim and S_rim) is within the 
uncertainty of the CXB surface brightness. We conclude 
that no additional power law component is required for 
our fit. 

3. 3. Model Fit for the Inner Regions ' Spectra 

We also applied the single-component VNEI model for 
each inner region. However, this model was inadequate 
to fit the data due to its simplicity. As previously shown 
by Katsuda et al. (2009), the thermal emission from the 
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Table 3. Spectral parameters 



1.56x10 



-11 



Component Parameter Value 

NW_riin E_rim 

CXB model 

absorption Nh [10^^ em~^] ... 3.6 (fixed) 

Broken Power-law Ti 2.0 (fixed) 

Ebroak [keV] 0.7 (fixed) 

T2 1.4 (fixed) 

surface brightness* 1.88x10"" 1.61x10"" 

LHB model 

APEC kTe [keV] 0.1 (fixed) 

Abundances 1 (fixed) 

surface brightness* 3.3x10"^^ (fixed) 

VNEI model 

absorption Nh [lO^i cm^^] . . . 2.3±0.1 2.0±0.1 

kTe [keV] 0.45±0.01 0.49±0.01 

log(r [s cm-3]) . . . 10.64±0.02 10.54^° °^ 

N 0.12±0.01 0.15±0.01 

O 0.24±0.01 0.22±0.01 

Ne 0.42±0.01 0.44±0.01 

Mg 0.37±0.02 0.34±0.01 

Si (=S) 0.71±0.08 0.32±0.01 

Fe (=Ni) 0.36±0.01 0.37±0.01 

xVdof 1395/1245 1441/1190 

* 2.0-10.0 keV surface brightness in units of erg cm"'^s"^deg"'^ 

Other elements are fixed to solar values. 

The errors are in the range Ax^ < 2.7 on one parameter. 



3.1±0.1 
0.41±0.01 

0.50±0.01 
0.25±0.01 
0.57±0.01 
0.47±0.01 
0.48±0.14 
0.40±0.01 

1183/1031 




Fig. 5. NW_rim spectrum fitted with the single-component 
VNEI model. 

interior of G156.2-I-5.7 consists of two components with 
different origins; a high-temperature ejecta component 
and a low-temperature ISM component. Thus we applied 
a two-component VNEI model for the spectra obtained 
from Center, East and NW. In this model, we fixed the 
ISM abundances to the average values based on the fit- 
ting results of rim observations (see table 3). The av- 
erage abundances of the ISM component are set as fol- 
lows: N=0.23, 0=0.24, Ne=0.45, Mg=0.39, Si(=S)=0.50, 
Fe(=Ni)=0.38. We fixed the other abundances to the solar 
values. We also set kTe to be 0.45 keV which is an average 
electron temperature for the rim regions of G156.2-I-5.7. 
The values of r and emission measure were allowed to vary. 



Meanwhile, in the ejecta component, the abundances of 
O, Ne, Mg, Si, S and Fe were free while we set the abun- 
dances of C and N equal to O, Ni equal to Fe. All other 
abundances were fixed to their solar values. The other 
free parameters were kTe, t, emission measure, and the 
column density iVn. 

We note that the temperature of the low-kT component 
may be higher than that of the rim because the line-of- 
sight ISM component includes ambient medium shocked 
earlier in the evolution of the remnant when the shock 
velocity was higher. In general, higher ISM temperatures 
will increase the line emission from species such as Si and 
S, which may affect the ejecta abundance measurement. 
Since it is impossible to measure the real average tem- 
perature of the ISM component, we examined how the 
metal abundances change by increasing the ISM temper- 
ature from 0.45 keV to 2.0 keV. As a result, we concluded 
that the ISM temperature hardly affect the ejecta abun- 
dances within the range of the error. We thus fixed the 
ISM temperature to be 0.45 keV for the inner regions' 
spectra. 

4. Discussion 

Ejecta Structure 

Figure 6 shows band images of O Vll, O Vlll, Ne IX, 
Mg XI, Si XIII and S XV using all the Suzaku data. In figure 
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Fig. 6. Band images of several elements. The white contours were derived from a soft X-ray image taken by ROSAT PSPC. 

6, both the Si and the S distributions appear to be cen- 
trally concentrated, while those of OVII and OVIII show 
a clear limb brightening, which suggests that the observed 
heavy elements originated from the ejecta. Furthermore, 
figure 6 suggests that the distribution centers of Si and S 
shift toward the northwest. 

To examine the radial distribution of the plasma, we 
divided the data into annular regions, as shown in fig- 
ure 1 right. As the annular center we used 04'^58™54^.9, 
51°44'36".l (J2000), the optical axis during the observa- 
tion of Center. We used an annular width of 2', consis- 
tent with the angular resolution of the XIS. In order to 
estimate the boundary of the ejecta-dominated region, we 
applied both the single-component VNEI model and the 
two-component VNEI model to the spectra obtained from 
the E_rim and NW_rim. We determined which model is 
preferred by using the F-test with a significance level of 
99%. We found that the second component is required 
in the NW_rim only within roughly ^85% of the shock 
radius (^52'). Figure 7 shows example spectra obtained 
from the annular regions 25' away from the annular cen- 
ter toward the cast (E-25; left) and the northwest (NW- 
25; right). The best-fit parameters are shown in table 
4. While both spectra require the two-component VNEI 
model, the abundances of the ejecta component are signifi- 
cantly different particularly for Si and S. We also note that 
there is no positive evidence that G156.2-I-5.7 is an 0-rich 
SNR as claimed by Hudaverdi et al. (2010). Their argu- 
ment is that the south and east side of their FoV are soft 
X-ray dominated (0.4-1.0 keV), suggesting the abundant 
presence of O. However, we consider that such excesses 



are mainly due to the contribution of the low-temperature 
ISM component. 




0.5 1 2 5 

Energy (keV) 




Energy (keV) 

Fig. 7. Example spectra of two annular regions in 
G156.2+5.7. Both were fitted with the two-component VNEI 
model. 

Figure 8 shows the radial profiles of the emission mea- 
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Table 4. Spectral parameters (see figure 7) 



Component 



Parameter 



Value 



-11 



Region E-25 Region NW-25 

CXB model 

absorption Nh [10^^ cm^^] 3.6 (fixed) 

Broken Power-law Ti 2.0 (fixed) 

Ebrcak [keV] 0.7 (fixed) 

r2 1.4 (fixed) 

surface brightness* 1.48x10^^^ 1.74x1 

LHB model 

APEC fcTe [keV] 0.1 (fixed) 

Abundances 1 (fixed) 

surface brightness* 3.3x10"^^ (fixed) 

VNEI model 

absorption Nr [10^^ cm^^j 1.8±0.1 

Low-kTf, component kT^ [keV] 0.45 

iog(T [s cm-3]) ii.ssiEj-;;? 

N 0.23 

O 0.24 

Ne 0.45 

Mg 0.39 

Si 0.50 

S 0.50 

Fe (=Ni) 0.38 

High-kTe component kT^ [keV] 0.59j:!^:[ji 

log(T [s cm-3]) 11.29±0.04 

O (=C=N) l-18io'2? 

Ne 0.51±0'.ll 

Mg 0.83i°:i^ 

Si 0.48±0.22 

S <1.5 

Fe (=Ni) 1.05±0.07 

xVdof 379/326 

* 2.0-10.0 keV surface brightness in units of erg cm~'^s~^deg^'^ 
Other elements are fixed to solar values. 
The errors are in the range Ax^ < 2.7 on one parameter. 



3.4±0.1 

(fixed) 

iu.uu_o.o8 

(fixed) 

(fixed) 

(fixed) 

(fixed) 

(fixed) 

(fixed) 

(fixed) 

0.58±0.01 
n 14+"°^ 

^^■-'^^-0.04 

0.61±0.09 
0.56±0.08 
0.63±0.08 
2.25±0.20 
3.40±0.96 
1.45±0.06 

533/453 



sure for various metals (O, Ne, Mg, Si, S, and Fe) in the 
ejecta. We found prominent central concentrations for Si, 
S, and Fe, while the lighter elements (O, Ne, and Mg) arc 
distributed more uniformly. The emission measure of Si 
peaks '^10' from the annular center and the overall dis- 
tribution is quite similar to that of S. The co-existence 
between Si and S is consistent with a theoretical picture 
that they are mainly produced in the static 0-burning 
layer during the stellar evolution and the incomplete Si- 
burning layer during the SN explosion. Similar shifts are 
observed in other evolved SNRs, for example, the Cygnus 
Loop (Uchida et al. 2009). If such inhomogeneities are 
real, one of the possible causes is an asymmetric explosion 
of a progenitor star. While many SN models have been 
proposed in recent decades, a model invoking delayed ex- 
plosion caused by neutrino heating is worth considering in 
the context of these data. For example, Marck & Janka 
(2009) simulated an SN explosion of 15M0 caused by the 



SASI-aided core-collapse (SASI; standing accretion shock 
instability). They expect that the heavy elements, such 
as Si and Fe, are ejected in one direction, while the lighter 
elements, such as O, Ne, and Mg, are ejected in all di- 
rections. However, since our FoV is limited, we cannot 
determine the precise distribution center of each ejected 
element. Also, the annular center does not necessarily 
correspond to the explosion center. It should be care- 
fully studied whether the origin of such inhomogeneous 
distribution is actually the result of an asymmetric SN 
explosion or not. 

^.2. Origin of hard X-ray emission 

As shown in section 3.2, we found no robust evidence 
for non-thermal emission in the NW_rim when we ex- 
cluded the region contaminated by the cluster 2XMM 
J045637.2-I-522411. The consistency of the measured 
CXB surface brightness from the NW_rim with other ob- 
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Fig. 8. Distributions of emission measure for various metals (O, Ne, Mg, Si, S, and Fe) in the ejecta. Open circle, cross, square, 
and triangle correspond to the spectral extracted regions. East, Center, NW, and NW_rim, respectively. The horizontal axis shows 
a distance from the center of Center. The positive and the negative values represent directions toward the northwest and the east, 
respectively. 
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servations' measurements (see table 2) as well as with val- 
ues found in the literature strengthens the conclusion that 
there is no non-thermal emission from the NW_rim. 

We next investigate spatial fluctuations of the hard X- 
ray surface brightness by using all the spectra obtained 
from the annular regions described above (see figure 2). 
From figure 2 left, there is no correlation between the 
radio synchrotron and the hard X-ray emission contrary 
to the prediction by Katsuda et al. (2009). The hard 
X-ray distribution does not have a bilaterally symmetric 
morphology like some non-thermal SNRs (e.g., SN1006; 
Koyama et al. 1995). Figure 2 right shows a histogram of 
the calculated surface brightnesses. The horizontal arrow 
represents the 1-a spatial fluctuation of the CXB surface 
brightness estimated by Tawa (2008). Surface bright- 
nesses in almost all regions (^94%) are within the 3-ct 
fluctuation, which suggests that the observed hard X-ray 
emission is dominated by the CXB. Although Katsuda 
et al. (2009) argued that non-thermal emission was also 
detected at the SNR's center, figure 2 shows that the hard 
surface brightness from this region is also consistent with 
CXB fluctuations, as is the case in the NW_rim. As a re- 
sult, there is no robust evidence for non-thermal emission 
from this remnant. 

4.3. Estimating the age of G156.2+5.7 

As shown above, the observed hard X-ray emission is 
dominated by the CXB everywhere except in the region of 
4'.0 around 2XMM J045637.2-H522411. This fact supports 
the idea that G156. 2-1-5.7 is an evolved SNR, despite the 
claim of Gcrardy & Fesen (2007) that it was a young 
nearby one. As explained in section 4.1, G156.2-I-5.7 has 
a layered metallicity structure that is observed in other 
evolved SNRs (e.g., Cygnus Loop; Tsunemi et al. 2007). 
Such structure suggests that the reverse shock has already 
reached at the center of G156. 2-1-5.7. Furthermore, the 
electron temperature (< 1 keV) is much lower than that 
in many young SNRs. All the results support the picture 
of an evolved SNR. For these reasons, we prefer ~1.1 kpc 
(Pfeffermann et al. 1991) to ^0.3 kpc (Gerardy & Fesen 
2007) for its distance d in the following analysis. 

Assuming thermal equilibrium between electrons and 
protons, we find that the forward shock velocity Vs is esti- 
mated to be Vs ~ 614(fcTe/0.45keV)°-^ km s~^. However, 
Ghavamian et al. (2001) calculate an initial equilibration 
of only 50% for shock speeds of ~ 600km s^^. In such a 
non-equilbrium case, the electron temperature kTg, and 
the proton temperature, kTp are expressed in the follow- 
ing equations (Ghavamian et al. 2001): 

1 3 

kTpC^-—[mefcct + {2-fcq)mp]v^^, (1) 
1 3 

fc^e ~ - — [TOp/cq + (2 - /cq)me] , (2) 



where m^^e are the proton and electron masses, respec- 
tively. If the initial electron-proton equilibration, f^q is 
0.5, the forward shock velocity is 960{kTe/0A5keY)°-^ 



km s^^, which gives an upper limit of Vg. 

On the other hand, the ambient densi- 
ties (no) for the NW_rim and the E_rim are 
n^o.NW = 0.084(rf/l.l kpc)~^ cm"-^ and no,E — 
0.091((i/l.l kpc)^^ cm~^, respectively, by fitting their 
emission measure profiles with a Sedov model (Katsuda 
et al. 2009). Due to the lack of statistics, it is hard to 
estimate the ambient density for the S_rim in the same 
way. We therefore consider the average value of two cal- 
culated no to be a global ambient density of G156.2+5.7: 
hence no = 0.088(c?/l.l kpc)~^ cm""^. Applying a simple 
Sedov analysis, the age t^ = t/W^ yr is as follows: 

' V614kms-i/ VlO^^ergy V0.088cm-3/ ^ ^ 
for the thermal equilibrium case, and 

^~ V960kms-i/ VlO^'ergy Vo.088cm-3y' 

for the non-equilibrium case. Here Eq is the explosion 
energy. As a result, we prefer the estimate from the 
ASCA observation (15000 yr; Yamauchi et al. 1999) to 
that from the ROSAT observation (26000 yr; Pfeffermann 
et al. 1991). In any cases, G156.2-f 5.7 is not a young SNR, 
as claimed by Gerardy & Fesen (2007). While Katsuda 
et al. (2009) pointed out that non-thermal X-ray emission 
would be detectable even from such an evolved SNR if 
the flux ratio of non-thermal to thermal emission is higher 
than those in a typical evolved SNR, it is more reasonable 
to conclude that G156.2-I-5.7 is, itself, a typical evolved 
SNR, since CXB exceeds that of any X-ray synchrotron 
component. 



5. Summary 

We performed a set of 6 pointing observations of 
G156.2-I-5.7 with Suzaku. The spectra were well fit- 
ted with two-component or single-component VNEI (with 
CXB and LHB) models for the inner regions and rim re- 
gions, respectively. Our results are consistent with the 
conclusion of Katsuda et al. (2009), that these compo- 
nents represent reverse shock-heated cjecta and forward 
shock-heated ISM, respectively. 

From the ejecta component, we found prominent central 
concentrations of Si, S, and Fe, while the lighter elements 
(O, Ne, and Mg) are distributed more uniformly. Such 
distributions reflect a layered-metallicity structure of the 
progenitor star. 

A single-component VNEI model provided an accept- 
able flt for all the rim regions; no additional power law 
component was required for our fit. We found no ro- 
bust evidence for the non-thermal emission in the NW_rim 
when we excluded the region contaminated by the cluster 
of galaxies, 2XMM J045637.2-h522411. Accordingly, we 
conclude that the hard X-ray emission in G156.2-I-5.7 is 
sufficiently explained as CXB emission. 

The estimated forward shock velocity is relatively 
slow compared with young SNRs: Vs ~ 614 — 
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960 (fcTe/0.45keV)''-5 km s"^ We estimated its age to 
be 7,000-15,000(d/l.lkpc)-i yr; therefore, G156.2+5.7 is 
not a young, but an evolved SNR. 
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